Charge-induced conformational changes of dendrimers 
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We study the effect of chargeable monomers on the conformation of dendrimers of low generation 
by computer simulations, employing bare Coulomb interactions. The presence of the latter 
leads to an increase in size of the dendrimer due to a combined effect of electrostatic repulsion 
and the presence of counterions within the dendrimer, and also enhances a shell-like structure 
for the monomers of different generations. In the resulting structures the bond-length between 
monomers, especially near the center, will increase to facilitate a more effective usage of space in 
the outer-regions of the dendrimer. 
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I. INTRODUCTION 

Dendrimers have been the subject of intensive inves- 
tigations ever since their synthesis 1 in the late 1970s. 
They are characterized by a high degree of monodisper- 
sity and a well-defined, highly branched internal struc- 
ture; efficient dendrimer assembly has been boosted by 
recent progress in synthetic techniques P A great deal of 
research activity has focused on the issue of whether they 
possess an open, dense-shell or a collapsed, dense-core 
configuration, the motivation arising by the potential to 
employ them as hollow, carrier-type molecules in the for- 
mer case. For neutral dendrimers, a large number of sim- 
ulation studies,^ careful self-consistent field calculations^ 
and not least scattering experiments^' have revealed that 
the dense-shell conformation is not the real one. Due 
to back-folding of the end-groups, caused by entropic 
considerations, a dense-core calculation results instead, 
leading even to compact, ha rd-s phere- like conformations 
at high generation numbers)^ From the point of view 
of applications, this may sound like a disappointing re- 
sult, as one would like to have dense-shell molecules. 
However, seen from the angle of fundamental research, 
the growing compactness of dendrimers with increasing 
generation number is very welcome, since it allows to 
use the m as model colloidal/nano particles with tunable 
stifmessP^H bridging the gap between flexible polymers 
and rigid spheres. 

The issue of dendrimer conformations is less clear when 
charged or polyelectrolyte dendrimers are considered. 
Charge on the building blocks of, e.g., poly(amidoaminc) 
(PAMAM) dendrimers can be manipulated by changing 
the pH of the solution^ and their conformations can be 
further influenced by added salt. The expectation that 
charged dendrimers may achieve stretching is based on 
experience with other branched polyelectrolytes (PE), 
such as PE-stars, for instance.^ However, the SANS- 
study of Nisato et alphas led to a negative result: the 
size of dendrimers is insensitive to pH changes and thus 
to charge. This experimental fact is at odds with the ear- 
lier work by Welch and Muthukumarp^ w ho predicted, 
by means of simulation and theory, an 'opening up' of PE 



dendrimers upon increase of the number of charged units; 
similar conclusions were reached in the Brownian Dy- 
namics simulations of Lyulin et However, both sim- 
ulational works quoted above employed a Debye-Hiickel 
(screened Coulomb) interaction potential acting between 
charged units, treating thereby the counterions as a con- 
tinuum. Giupponi et al., on the other hand, pointed out 
that a more realistic treatment of the molecules should 
employ the bare Coulomb interactions among all involved 
species (monomers and counterions) , which is the method 
they employed in their own, recent simulations.^ Inter- 
estingly enough, they found that the dendrimer size is 
indeed very weakly dependent on charge, due to local 
charge neutrality, a condition that is masked when one 
employs the Debye-Hiickel approximation. 



Though the work of Giupponi et alf^ has illuminated 
a number of issues pertaining to monovalently charged 
monomer units accompanied by monovalent counteri- 
ons, the question of the influence of valency on den- 
drimer conformations has not been studied so far. The 
purpose of this work is to examine precisely this issue, 
which appears relevant on the grounds that counterion 
valency is known to bring about drastic changes in, e.g., 
the size of spherical PE brushes!^ We focus thereby 
on dendrimers of the fourth generation (G = 4) and 
examine separately the conformations of neutral den- 
drimers, as a reference point, as well as the combina- 
tions (Z m ,Z c ) = (1,1), (1,2), (2,1) and (2,2), where Z m 
stands for the valency of the monomers and Z c for the 
valency of the counterions. We find that the case of diva- 
lent monomers and monovalent counterions brings about 
a substantial change of the dendrimer size, accompanied 
by a strong stretching of the chemical bonds, whereas the 
same phenomena are less pronounced for the other two 
cases. After describing our simulation model in Sec. II, 
we present and discuss the results in Sec. Ill, whereas in 
Sec. IV we draw our conclusions. 
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II. THE SIMULATION MODEL 

The dendrimers wc use within our simulations are built 
from a central pair of joined monomers, the so-called gen- 
eration 0. A successive generation g + 1 of dendrimer is 
formed by connecting two additional monomers, the func- 
tionality of the dendrimer is therefore three, to each outer 
monomer of the dendrimer of generation g. In doing so, 
the number of monomers n(g) of a given generation g in 
a dendrimer follows a simple power law, i.e. n(g) = 2 9+1 . 

There arc in general three types of interactions be- 
tween monomers. The first type of interaction prevents 
the collapse of monomers onto each other and is a short- 
range repulsive interaction given by a simple, shifted and 
purely repulsive Lenard-Jones potential 




r <r c 
r > r r 



(1) 



where a and e are the unity of length and energy respec- 
tively, and r c = 2 1 / 6 er is the range of the interaction. 

The second type of interaction is of an attractive na- 
ture and describes the bonds between joined monomers 
in order to prevent the molecule from flying apart. This 
interaction is described by a FENE potential 
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where Rq = 1. 5<r is the maximum allowed distance be- 
tween two connected monomers. 

The last type of interaction to be included in the model 
is the Coulomb potential between charged monomers 



Vbouiomb — ksTXs 1 J 
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with Zi and Zj the charge numbers, &b the Boltzmann 
constant, T the temperature, and Ab the Bjerrum length 
given by 
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We have performed molecular dynamics simulations at 
constant d ensity and temperature using the Nose-Hoover 
thermostat! 18 * 19 ! on dendrimers of generation 4, i.e., a 
dendrimer formed by 62 identical monomers, where the 
monomers with charge number Z m are either neutral, 
monovalent, or divalent. In order to guarantee charge 
neutrality the charged dendrimers need to be balanced, 
for which we in either case used both monovalent and 
divalent counterions (charge number Z c ). Apart from 
the Coulomb interaction between monomers and counte- 
rions, as well as between the counterions themselves, we 
use for simplicity the same short range repulsion ([lj as 
is used for the monomers. 

Simulation parameters are chosen such that our unit 
of length a = 2.84A and we fixed the temperature to 



T = 1.2e/fce. Using an implicit solvent that mimics 
the behavior of water, i.e., a temperature of 300K and 
a relative permittivity e r = 80 we arrive at Ab/ct = 3. 
Periodic boundary conditions have been applied in com- 
bination with the Ewald summation method to include 
the long-range electrostatic interactions. The volume of 
the simulation box was chosen such that effectively the 
dendrimers can be considered to be independent, indi- 
cated by the independence of the results on larger vol- 
ume sizes. The initial configurations were equilibrated 
over times long enough for the counterions to diffuse in 
to the core of the dendrimer, and to reach a steady state 
for a in- and out-flux. The simulation time is chosen long 
enough for the individual counter-ions to explore the full 
dendrimer, i.e., inner and outer ranges. 



III. RESULTS 

In order to characterize the size of the dendrimer, the 
radius of gyration R g is measured, which is defined by 
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where the summation runs over the N positions f\ of 
the monomers of the dendrimer and r^u is the center of 
mass of the dendrimer. 

The resulting values for the neutral and charged den- 
drimers are listed in Table HI The first observation we 
can make, is that by charging the dendrimer the radius 
of gyration increases. This is actually not so surpris- 
ing, because the monomers will repel each other due to 
their charges. What is more interesting is the fact that 
the behavior on charging is not monotonic, i.e., the ra- 
dius of gyration of the divalently charged monomers and 
counterions is smaller with respect to that of the mono- 
valent case. This suggests that the swelling of the den- 
drimer due to the increasing charge is counteracted by 
the transport of counterions into the deeper regions of 
the dendrimer, w hich re sults in screening effects and re- 
duces the swelling. ^ 1 21 ' 22 1 



Z m — Zc 
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R g /a 


neutral 


6.31 ± 0.02 


2.51 ± 0.01 


1-1 


8.20 ± 0.02 


2.86 ± 0.01 


1-2 


7.12 ± 0.02 


2.67 ± 0.01 


2-1 


9.27 ± 0.02 


3.04 ± 0.01 


2-2 


7.49 ± 0.02 


2.74 ± 0.01 



TABLE I: The radius of gyration for the neutral dendrimer 
and the dendrimers with charged monomers and counterions, 
with charge numbers Z m and Z c respectively. 

It also appears that divalent counterions reduce the 
size of the dendrimer with respect to that of monovalent 
counterions. Not only do they result in a better screen- 
ing of the monomers, but also less of them are required 
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which also leads to less steric hindrance. Using a similar 
argument, the monovalent monomers reduce the size of 
the dendrimer with respect of that of divalent monomers, 
because in the latter case more counterions are required 
within the dendrimer, preventing its shrinking. Note 
that for a neutral dendrimer there is no counterpart for 
the counterions present within the dendrimer, hence it 
is only the short-range repulsion between the monomers 
that prevents the dendrimer from collapsing. 

The radius of gyration is also experimentally observ- 
able from the small wave-vector limit of the form factor 
F(q) defined by 




FIG. 2: Radial density profiles for the monomers (a) and 
counterions (b), measured with respect to the center-of-mass 
of the dendrimer. The curves in the left-hand figure have been 
shifted to facilitate a comparison. 
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By expanding in small wave-vectors q and averaging over 
the directions it is easy to show that this results in 



F(q) = N 
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In the present case this expansion is valid up to qa ~ 0.2 
as can be seen in Fig. [I] where we compare the measured 
form factor, Eq. ([6|, from the simulation with the small 
wave- vector limit, Eq. Q. 
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FIG. 1: Comparison of the measured form factor, Eq. Q 
from the simulation (sim), with the small wave- vector limit, 
Eq. Q, denoted (th). 

In Fig. [2] the radial density profiles of the monomers 
and counterions are shown, both measured with respect 
to the center of mass of the dendrimer. Note that at small 
distances from the center of mass the noisy behavior is 
purely due a poor statistical sampling caused by the lack 
of particles in that region. 

It is immediately clear that the presence of charge on 
the monomers results in a much more structured den- 
sity profile, in which the monomers are mostly found 
in a shell- like structure. The mutual repulsion of the 
monomers due to their charge not only leads to a larger 
size of the dendrimer but it enables the counterions to 



diffuse into the dendrimer as well; this can be seen from 
their density profiles. Their distribution, however, shows 
less structure and has a wider range. The latter can also 
easily be understood, since the presence of counterions 
within the dendrimer will also lead to steric hindrance 
and there is an obvious entropic gain in surrounding the 
dendrimer and move in the region lying outside its extent. 




FIG. 3: Radial density profiles of the monomers and coun- 
terions for the monovalently [(a),(b)] and divalently [(c), (d)] 
charged monomers, and monovalently [(a),(c)] and divalently 
[(b), (d)] counterions. 

To clarify the ordering in the charged dendrimer cases, 
the same data are shown in Fig. [3] where the radial den- 
sity profiles of the monomers and counterions are directly 
compared. Although the highest density for the counte- 
rions is reached roughly in between the first and second 
peak of the monomer distributions (except for the case 
Z m = 1 and Z c = 2) , there is no layered structure present 
for the counterions. This is also confirmed by the absence 
of plateaus in the cumulative counterion density profile 
(not shown). 

The absence of charges on the monomers in the neu- 
tral dendrimer does not only affect its size, i.e., it is more 
compact, but it also modifies its internal structure as is 
illustrated in Fig. [4] where the density profiles for the 
different generations are plotted versus the distance to 
the center of mass. This reveals that the density of gen- 
eration 4 monomers at smaller distances is larger than 
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that of those from generation 3. In other words, the den- 
drimer starts to fold in to itself. Whereas the repulsive 
Coulomb interactions between charged monomers tend to 
stretch the dendrimer causing a loss in entropy, the neu- 
tral monomers just feel the short-range Lennard- Jones 
repulsion and can exploit the open space in the core of 
the dendrimeric structure. A not unimportant other rea- 
son is the absence of the equivalent of counterions within 
the dendrimer. This effect is only weakly visible in the 
cases of the divalently charged counterions, but it is to be 
expected that it will be more pronounced for increasing 
value of the generation of the dendrimer. 
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FIG. 4: Radial density profiles of the monomers in the neutral 
dendrimer decomposed in the contributions stemming from 
each generation. 

A more detailed description of the internal structure 
of the dendrimer can be obtained by analyzing the the 
bond-lengths b, shown in Figs. [5] and [6] The bond- length 
probability distributions, P(b), are decomposed per gen- 
eration, whereby a bond of generation n is formed by a 
monomer of generation g with its parent of generation 
g — 1 (with the exception of generation 0). For the neu- 
tral dendrimer only a small shift in the distribution of 
lengths is found towards shorter bond-lengths for higher 
generations, i.e., bonds near the center of the dendrimer 
tend to be more stretched than those near the border. 
This implies that there is a collective behavior in which 
the mutual repulsion of monomers in higher generations 
that prefers to expand the dendrimer, forcing the fewer 
central bonds to stretch. 

Fig. [6] shows the same distributions but now for the 
charged dendrimer cases. The first thing one can ob- 
serve is that the bond-lengths in the charged cases are 
more stretched. This is not so surprising, since the bond- 
length is directly affected by the mutual charge repul- 
sion of the monomers and even more so for the divalent 
monomers. Also the tendency for the central bonds to 
be more stretched than those in the outer regions is ap- 
parent. The most interesting case, however, is that of 
the divalent monomers and monovalent counterions, in 
which case the stretching of the bond-lengths shifts sig- 
nificantly with the generation. In Fig.[7]this is illustrated 
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FIG. 5: Probability distribution of bond-lengths for each gen- 
eration for the neutral dendrimer. 

even more clearly by comparing the bonds of generation 
and 4 for the various models. This suggests that the 
abundance of monovalent counterions required for an ef- 
fective screening of the divalent monomer bonds can not 
be packed within the core of the dendrimer. 
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FIG. 6: Probability distribution of bond-lengths for each gen- 
eration for the for the monovalently [(a),(b)] and divalently 
[(c), (d)] charged monomers, and monovalently [(a),(c)] and 
divalently [(b), (d)] charged counterions. 




FIG. 7: Probability distribution of bond-lengths of generation 
(a) and generation 4 (b) for the different cases monomers 
and counterions. 

In the final method we used to examine the internal 
structure of the dendrimer, we consider the angles be- 
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FIG. 8: Schematic representation of the local structure of a 
monomer inside a dendrimer indicating the different 9 and 4> 
angles mentioned in the text. 



tween bonds in the dendrimer. The angle <f> as illustrated 
in Fig. [8] is the angle between the bonds of a monomer 
of generation g and the two monomers bounded to it of 
generation g + 1 (a bond generation g + 1). The angle 8 
is the angle between a bond of generation g and a bond 
of generation g + 1. 




Neutral 

- z m =\,z c =\ 

- — Z* = l,Z c = 2 

- — Z™ = 2,z'=l 
— Z m = l,Z t = 2 








/ / w 




/ \w 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
$1% 



FIG. 9: Probability distribution of the angles 9 [(a)] and <j) 
[(b)] for the different dendrimer cases averaged over all nodes 
in the dendrimer. 

The results of the angle analysis are shown in Fig. [9] 
where we have averaged over all 30 nodes in the den- 
drimer. In the case of the neutral dendrimer the prob- 
ability distributions are broader than in the case of the 
charged dendrimers. This is partially due to the fact 
that there is no penalty for smaller angles, whereas in the 
charged dendrimer cases the angles will be more focused 
in the equal angles area due to the Coulomb repulsion. 




FIG. 10: Probability of the angle distributions of a dendrimer 
with monovalent monomers and counterions. (a): angle 9; 
(b): angle (j>. 



Also these angle distributions will, just like in the 
case of the bond-lengths, depend on the position of the 
monomer within the dendrimer, i.e. its generation, as can 
be seen in Fig. [10] in the case of monovalent monomers 
and dendrimers. For nodes near the core the <fi angles are 
more or less centered about 120°, but for nodes near the 
boundary of the dendrimer the peaks shifts to smaller 
angles. This is specially clear for the angle <f> for the 
outer most nodes lies between 70° and 80° for the vari- 
ous cases presented here. At the same time the 8 angles 
show a similar behavior, except that they become larger, 
albeit that the shift is not as significant as is the case for 
0. 

This is not the only information that can be extracted 
from the angles 8 and (p. In general the three bonds at a 
given node will not lie within a single plane, but span a 
space angle, also called spherical excess,^ e, which ranges 
from zero to 2ir. Both extremes can only be obtained in 
the limiting case of three bonds lying within the same 
plane. The fact that for such a planar configuration the 
spherical excess can take two values can be understood 
by realizing that two cases need to be distinguished. If 
the sum of the three angles between the pairs of bonds 
9i + 8 2 + 4> = 2/T, the space angle reaches it maximum 2tt 
and spans a half sphere. In the other case the largest of 
the three bond angles equals the sum of the smaller ones, 
e.g. 8i + 82 — <fi and the spherical excess reaches it mini- 
mum 0. In other words, when the three bonds completely 
open up and stretch away from one another, one has the 
maximum value, e = 2tt; on the other hand, when there 
is complete backfolding of the two bond connecting the g- 
monomer with the two g+ 1-monomers towards the bond 
between the g-monomcr and the g — 1-monomcr in Fig. 8, 
then e = 0. In this fashion, the probability distribution 
P(e) of the spherical excess offers valuable information on 
the presence of back-folding within the dendrimer. A dis- 
tribution with strongly suppressed values around e = 
and high values around e = 2n implies a stretched den- 
drimer with planar three-bond junctions at the branching 
points. Relatively flat distributions with non-negligible 
values around e — point rather to strongly back-folded 
dendrimers, akin to neutral ones. 

The probability distribution of the spherical excess 
for the neutral dendrimer and for the case of a charged 
dendrimer with both monomers and counterions mono- 
valcntly charged, is shown in Fig. [TT] It reveals that 
the spherical excess decreases for nodes that lie further 
away from the core of the dendrimer. The preference 
for spherical excess to be large near the center of the 
dendrimer can be understood by realizing that such be- 
havior will maximize stretching in the dendrimer and en- 
able to exploit the so available space for the more distant 
monomers more effectively. The decrease of the spherical 
excess is particular striking in the case of the neutral den- 
drimer and is a necessary consequence from our earlier 
observation that the monomers from the last generation 
can fold back into the dendrimer.^ As can be seen in 
Fig. |ll[ b), however, and in conjunction with Fig. 11 'a), 
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the Z m = Z c — 1-dendrimer shows a considerably higher 
degree of stretching than the neutral one; its distribution 
of the spherical excess is shifted towards higher values 
and only at the third generation do we observe a nonva- 
nishing contribution around e = 0. 

Fig. [12] shows the spherical excess of the nodes of gen- 
eration and 3 for the various models discussed. In 
full agreement with all our previous measures, we find 
that charged dendrimers carrying divalent monomers and 
monovalent counterions show the most dramatic degree 
of stretching, all the way from the inner to the out- 
ermost generations. Physically, this is caused by the 
strong attraction between the monomers and the counte- 
rions. In this case, there are twice as many counterions 
as in the cases Z m = Z c = 1 and Z m = Z c = 2 that 
have to be accommodated in the inner region of the den- 
drimer, giving rise to stretching to create space. Two 
counterions per monomer are required to achieve local 
charge neutrality,^ rendering the composite, monomer- 
countcrion groups too big to allow for the backfolding 
seen for neutral dendrimers. 
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FIG. 11: Spherical excess for the neutral dendrimer [(a)] and 
monovalently charged monomers and counterions [(b)]. 




chargeable monomers will have on their size and internal 
structure. In particular, we have allowed for the possi- 
bility of divalent monomer units and/or counterions and 
monitored the conformational changes induced by the 
various scenarios. For the simplest case, Z m = Z c = 1, 
our findings are in full agreement with the recent ones of 
Guipponi et al.\^ in the sense that no dramatic changes 
in the dendrimer conformation and size have been de- 
tected. This is also in agreement with the experimental 
results of Nisato et alf^ and at odds with the prediction 
of Welch and MuthukumarP^ At this point, it appears 
that the assumption of linear counterion screening, which 
has been employed in Rcf. [TS] in the form of a Debye- 
Hiickel approximation, is not valid when strong coun- 
terion absorption and condensation effects are present. 
More so at the nanoscale-distances involved in the inner 
of a dendrimer. 

At the same time, we could establish that a remark- 
able change of the size of the dendrimer, accompanied by 
a stiffening and stretching of its bonds, takes place when 
the former carries divalent chargeable groups, which re- 
lease two monovalent counterions per site. Here, al- 
though the resulting configuration is not yet of the dense- 
shell type, the dendrimer size grows by almost 50% in 
comparison to a neutral dendrimer and strong correla- 
tions between the monomer- and counterion profiles ap- 
pear. It is now reasonable to assume that, given a guest 
molecule that carries the same type of charge as the coun- 
terions and being small enough to 'fit' inside the stretched 
dendrimer, absorption will take place, as in this fashion 
a number of counterions will be released from the inner 
part of th e dendrimer, r esulting in a concomitant en- 
tropy gain! 24 | 25 | 26 | 27 | 28 | 2 S l Another possibility is to endow 
the dendrimers with a spacer length s > 1, which allows 
for creation for more space in their interior but at the 
same time it also increases the possibilities to backfold- 
ing. Work along these lines is currently in progress. 
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FIG. 12: Spherical excess of generation [(a)] and generation 
3 [(b)] nodes for the different cases monomers and counterions. 



IV. DISCUSSION 
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